A K 



4-k 



A N 



4_U 



^ 

^"7 ^"7 ^"7 \~7 \~7 \~7 



A N 



# 



4 



4 



/I K 



N 
4 



A 



A 



V N 



4 



b. 



7 



7 



7 



\'7 



A 



4 



4 



v 



7 



7 



A N 



^~7 




4 



4 



4 



4 



K 4 



7 ^ 



4 



V 



4 



k A 



7 ^ 



A 



A 



A N 



17 



/I N 



4 



V 



V ^ 



4 



4 



A 



A 



P 7 ^ 



\ 



N 



K A 



7 ^ 



K A 



V ^ 



K 4 



4 



4 



A 



A 



A 



"7 V 



7 



4 



4 



A 



v N 



4 



A 



A 



\ 



bw 



A 



7 ^ 



N 
4 



k A 



V N 



17 



N 
A 



7 



K A 



7 ^ 



7 



7 



7 



A N 



7 



/I K 



A 



k A 



V N 



>1 



4- 




K 



4 



7 ^ 



k A 



4 



k 4 



7 ^| 



k A 



7 ^ 



7 



V 



A 



V 

36. JOHN 



A K 
A U 



MANI 
COMflfrERTECHNJ 



4 



4 



7 



7 



4 



7 



4 



4 




DUE 
CUM! 



HOTES 
M.OGY DEPT 



7 



4 



7 ^ 



4 



4 



7 N 



^ 4 



N17 



K 4 



7 ^ 



4 



K 4 



7 ^ 



4 



7 



7 



7 



bw 



7 



4_K 4-bw 4_J^ 

*=# ^ 



7 



4-t^ 



N 4 



7 N 



4 



k 4 



4 



7 ^ 



7 



K 4 



7 ^ 



4 



N 4 



A N 



7 



K 4 



7 



4 



|7 



4 



b. 



7 



k 4 



7 



^"7 



4 



4 



N 4 



7 



7 



k 4 
K 4 



K 



7 



7 



*==t ^ 



N 

4 



k 4 



7 ^ 



4 



k A 



N 4 



V 



b, 



4 



4 



4 



7 ^ 



4 



7 ^ 



N 4 



7 n 



7 



K 4 



7 ^ 



7 



7 



/I K 



7 



K 4 



7 N 



4 



7 



A 



7 ^ 



4 



A. 

'7 ^ _ 



7 



V 



4 



4 



N 4 



7 ^ 



4 



k 4 



4 



7 



4-K 



V 



K 4 



7 .'SJ 



bw 



7 



4 



4 



7 ^ 



7 



K 4 



7 ^ 



b^ 



7 



h 4 



7 ^ 



7 



7 



7 



"7 



7 



NT 



4-K 



^"7 



4_K 



N~17 



b. 



7 



w 

A N 



7 ^v" 



7 



^"7 



"7 



DISCLAIMER 



Although each program has been tested by its contributor, no 
warranty, express or implied, is made by the contributor or 
1620 USERS Group, as to the accuracy and functioning of the 
program and related program material, nor shall the fact of 
distribution constitute any such warranty, and no responsibility 
is assumed by the contributor or 1620 USERS Group, in con- 
nection therewith. 



1620 USERS GROUP PROGRAM REVIEW AND EVALUATION 

(fill out iii typewriter or pencil, do not use ink) 

Program No. Date 

Program Name: 

1. Does the abstract adequately describe what the program is and what Yes No_ 

it does? 

Comment 

2. Does the program do what the abstract says ? Yes No_ 

Commen t 

3. Is the Description clear, understandable, and adequate? Yes No_ 

Comment 

4. Are the Operating Instructions understandable and in sufficient detail? Yes No_ 

Comment 

Are the Sense Switch options adequately described (if applicable)? Yes No_ 

Are the mnemonic labels identified or sufficiently understandable? Yes No_ 

Comment 

5. Does the source program compile satisfactorily (if applicable)? Yes No_ 

Comment 

6. Does the object program run satisfactorily? Yes No_ 

Comment 

7. Number of test cases run . Are any restrictions as to data, 

size, range, etc. covered adequately in description? Yes No_ 

Comment 

8. Does the Program Meet the minimal standards of the 1620 Users 

Group? Yes No_ 

C omment 

9. Were all necessary parts of the program received? Yes No_ 

Comment . 

10. Please list on the back any suggestions to improve the usefulness of the program. 
These will be passed onto the author for his consideration. 

Please return to: Your Name 



Mr. Richard L. Pratt Company 
Data Corporation 

7500 Old Xenia Pike Address 
Dayton, Ohio 45432 



User Group Code 



THIS REVIEW FORM IS PART OF THE 1620 USER GROUP ORGANIZATION'S PROGRAM 
REVIEW AND EVALUATION PROCEDURE. NONMEMBERS ARE CORDIALLY INVI TED 
TO PARTICIPATE IN THIS EVALUATION. 



11/09/fi4 



CROUT REDUCTION 



GROUT REDUCTION 



DECK KEY 



S. S. Millwright 



Deck 1 



Fortran Source deck 



Deck 2 



Object deck 



Deck 3 



Input to sample probli 



Modifications or revisions to this program, as they occur, 
will be announced in the appropriate Catalog of Programs 
for IBM Data Processing Systems. When such an announce- 
ment occurs, users should order a complete new program 
from the Program Information Department. 



Direct Inquiries to: 



Mr. R. C. Wenrick 

336 Woodward Road, S. E. 

Albuquerque, New Mexico 



1620 USERS Group Library 



Program Abstract 

Title: Crout Reduction 

Author; Organization: S. S. Millwright 

ACF Industries, Incorporated 
Albuquerque, New Mexico 

Direct Inquiries to: Mr. R. C. Wenrick Telephone: 247-0361, Ext. 507 

336 Woodward Road, S. E. 
Albuquerque, New Mexico 

Purpose/Description: Will evaluate determinants of maximum order 36, solve simultaneous 
linear systems with an indefinite number of constant vectors, and/or determine the inverse 
of maximum order 36 matrices. Singular input matrices are detected. Output includes 
data which indicates accuracy of solution vector. 

Mathematical Method: Crout Reduction. See Introduction to Numerical Analysis by 
F. B. Hildebrand, McGraw-Hill Book Co., Inc., 1956, pp. 429-435 

Restrictions, Range: Element in first row and first column of input matrix must not be zero. 

Storage Requirements: All of core. 

Equipment Specifications: 

Memory 20K 40K 60K_X_ K Automatic Divide: Yes No 

Indirect Addressing: Yes No Other Special Features Required_ 

Additional Remarks: This program was written in Fortran language. It was compiled on a 
machine equipped with the automatic divide/automatic floating point features. Smaller 
arrays could be specified in order to adapt this program to smaller machines. Successful 
solutions for systems up to ninth order have been accomplished with this program. 
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INTRODUCTION 

M"201 is o Fortran language program which employs the Crout Reduction ^ technique in 
order to evaluate determinants of maximum order 36, to solve simultaneous linear systems 
consisting of up to 36 equations in 36 unknowns with an indefinite number of constant 
vectors, and/or to calculate the inverses of maximum order 36 matrices. 

The coefficient matrix is tested for singularity by the program. If it is found to be singular, 
an error message is typed and punched and the program proceeds to the next problem with 
no programmed stop. 

If one or more solution vectors for simultaneous systems are desired, a matrix multiplication 
of the coefficient matrix times the solution vector will be performed to obtain a calculated 
constant vector. The difference between rhe input constant vector and the calculated 
constant vector for each solution vector will appear in the output. This difference will 
give an indication of the inaccuracy in the solution vector due to truncation. 

If the solution vector accuracy is found to be unsatisfactory, the user may find it desir- 
able to input the above defined "difference" vector as a new constant vector together 
with the original coefficient matrix. The new solution vector thus obtained may then 
be added to the original solution vector in order to obtain a better approximation of the 
true solution vector. This process may be repeated any number of times. The program 
does not perform this correction automatically; the input must be presented to the program 
as a separate problem and the "correction" vector added externally if this iteration 
technique is to be applied. 

The user must insure that the element in row 1, column 1 of the coeffltient matrix has 
a value other than zero. 

GENERAL ANALYSIS 

An augmented matrix of the original system of the form: 

1 

0) a,, a 12 ... a ln , c, 

a 21 a 22 ... a 2n J c 2 

a . a « . . . a„„ I c 

nl n2 nn 1 n 



Hildebrand, F. B.: "Introduction to Numerical Analysis", McGraw - Hill Book 
Co., Inc., New York (1956), pp429 - 435. 
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Where the a's are the coefficients of the unknowns and the c's are the constants can be 
reduced to an augmented matrix of the form: 



(2) 


a 'll a' 12 . . 


• a 'ln | c l 




a 21 a 22 • ' 


• a 2n j c 2 




a', a', . • 

nl n2 


a' 1 c' 
nn | n 



By application of the relations 



j-i 

alj= a y - I a; K a' Kj (i 



r 

aij-Ci/a-;-,) L aij - 1^ a^aVj 



First to find the elements of the first column (a'^ to a n1 ), next to find the remaining 
elements of the first row (a' to c"), next the remaining elements of the second column 
(a 22 to a^), and so on unfil the augmented matrix has been completely determined. 

It can be shown that the value of the determinant of the original coefficient matrix is 
equal to the product of the n number of elements of the reduced matrix which lie on 
the main diagonal . If I A \ denotes the value of this determinate then, symbol ical ly: 

IM= It a' u 
i = i 
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The solution vector for the original system can be obtained from the relation: 

n 

X; = c- - I a' lK x K 



The solution vector is calculated foot to head i.e., in the order X^j, . . . Xj 

It can also be shown that if the constant vector in (1) be replaced with an n* n order unit 
matrix and n solution vectors be calculated considering the n columns of the unit matrix 
to be n constant vectors proceeding from left to right, that the resulting n solution vectors 
are identically the left-to-right columnwise elements of the inverse of the coefficient 
matrix from (1). 



INPUT 

All input is from punched cards and shall consist of the following: 

Card No. Data CC Remarks 

1 Heading 1-55 Up to 55 Hollerith Characters 

2 N 1-5 15 - Order of matrix 
MTRX 6-10 15 - >O f invert matrix; 

^. O, do not invert matrix 
LSLN 11-15 15 - >O f solutions) to 

simultaneous system desired; 
^■O, no solution desired 
KNO 16-20 15 - number of constant vectors 

in input 

JVAL 21-25 15 -> O, evaluate deter- 

minant of coeff. matrix; ^ O, 
do not evaluate determinant. 

INXT 26-30 15- >0, read data for next 

problem; <0, exit from M-201 

3 to 2 + EN A(l, J) 1-60 6E10.4 - Values of coefficients 
(where E = of input matrix arranged by rows 
nearest integer > JN ") up to six elements per caraT 

k (elements must all be from the 



same row on any given card) 
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Card No. Data > CC Remarks 

3+ EN to C(l) 1-60 6E10.4- values of constant 

2 + E(N + KNO) vector elements up to six 

(where E = . per card (Elements must all 

nearest integer > N ) be from same constant vector 

k on any given card) 



OUTPUT 



All output is on punched cards except for the singular input matrix error message which 
will be typed as well as punched. 

The output deck is designed to be listed on the IBM 1403 printer using a Fortran output 
I ister program which skips to a new page for each new problem and deletes the Sequence 
numbers in cc 78 - 80 of the output deck. 

The output deck will consist of appropriate header cards, and any or all of the following 
data as called for on input card 2: 

1 . The value of the determinant of the coefficient matrix complete with an identifying 
label. 

2. The solution vector(s) together with the "difference" vector(s) which is the difference 
between the input constant vector and a calculated constant vector. These will 
appear with identifying header cards and labels. 

3. The inverse of the coefficient matrix listed by columns in five-column blocks. The 
column number will appear above the appropriate column; the raw number will appear 
to the left of the appropriate row. 



OPERATING INSTRUCTIONS 



1 . Load M-201 object deck followed by input deck into the read hopper of the IBM 
1622. 

2. Ready the punch hopper with blank cards. 



- 6- 



Set console switches: PROGRAM LISTING 

I/O to STOP 
Parity to STOP 

OFLO to PROGRAM The following pages contain the source language listing for M-201 . 

Sense switches — none tested. 

Depress RESET and LOAD. 
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C H-201 CROUT REDUCTION 

C HAY BE USED TO — 

C 1. EVALUATE DETERMINANTS {MAXIMUM 0R0ER=36) 

C 2. SOLVE UP TO 36 SIMULTANEOUS LINEAR EQUATIONS 

C 3. FIND THE INVERSE OF A MATRIX 

C fVE25? T IN FIRST R0W AN0 FIRST COLUMN, All, I), MAY NOT BE ZERO 

on 2.12! SI ° N M * »« 0, » AA<36 »36>'X<36),Ct36),CCC36),NQ(40> 
™9 PUNCH 1 I 

READ 3 
PUNCH 3 

READ l,N,MTRX,LSLN,KNO, JVAL, INXT 
DO 100 1=1, N 
DO 100 J=1,N,6 

.„« *"^ 2 »* (I » J, ' AlI »J + l>»*<I»J*2),A(I,J*3),A(I,J+4),A(I,J + 5) 

lOO CONTINUE 

00 200 1 = 1, N 
200 AA(I,l)=AtI,l) 

00 300 J=2,N 
300 AAC1, J)«A(l, J)/AU,1) 

00 400 1=2, N 

DO 400 J=2,N | ,, 

400 AA(I,J)=Q. / y: rj ,jt y u, .^^ifijy fQ 

J«2 
490 II-J 

DO "510 I«H,N 

LIMI*J-1 

00 500 K*1,LIM1 
500 AA(I,J)=AAtl,J)+(AA(I,K)«AA(K,J)) 
510 AAt I,J)=AU,J>-AA(I,J) 

IF(AACJ.J) 1520,900, 520 
520 IFtN-J 1700,700,530 
530 I*J 

J»J*1 

JJ=J 

00 610 J=JJ,N 
- LIK2-I-1 
00 600 K*1,LIM2 
600 «A<I,J)*«Af I,J)*tAA(I,K)«AA(K,J)) 
AAt I r J)*AII,J)-AA|I,J| 



610 AAtI,J)=AA(I,J)/AA(I,I) 

J*I + i 

GO TC 490 
700 IF( JVAL)800,800,720 
720 VALUE=AA(1,1) 

DO 710 1=2, N 
710 VALUE=VALUE«AAU, I) 

PUNCH 4, VALUE 
600 IF(LSLN)730,730,220 
220 00 230 1=1, N, 6 

READ 2,C(I),C(I+l),C(I+2),C(I+3),C(I+4),C(I+5) 
230 CONTINUE 

DO 240 1=1, N 

X(1 1=0. 
240 CC(I)=0. 

CCtl)=C(l)/AA(l,l) 

00 250 1=2, N 

LIM6=I-1 

00 260 K=1,LIM6 
260 CC(I)=CC(I)+(AA{ I,K)»CC(K)) 

CC( I l=C( I)-CC< I ) 
250 CC(I)=CC(I)/AA(I,I) 

X(N)=CCIN) 

LIM7»N-l 

00 270 I=1,LIM7 

II-N-I 

LIM8=II+1 

00 280 K=LIM8,N 
280 X(II)=X{II)+(AA(II,K)*X(K>) 
270 X(II)=CC(II)-X(II) 

00 290 1=1, N 

290 ccm=o. 

DO 291 1=1, N 
DO 291 J*1,N 

291 CC(II=CC(I)*(A{ I,JJ«Xl J) ) 
00 292 1 = 1, N 

292 CCtI)=C(I)-CC{I) 
PUNCH 8 

DO 293 1=1,40 



293 NG C I )=I 

CO 294 1=1, N 

punch 9,Ncm.xm,ccm 

294 CONTINUE 
KNO=KNC-l 

IFCKNC)730,730,220 
730 IFIKTRX)21C,210, 810 
810 C(1>=1. 

CC{1)=1./AA(1,1) 

DO 820 1=2, N 

C(I )=0. 

820 ccm = o. 

NIJ=N+4 

DO 830 1=1, NIJ 

CO 830 J=1,NIJ 
830 A(I,J)=0. 

J=0 
840 J=J + 1 

CO 860 1=2, N 

LIM3=I-1 

CO 890 K=1,LIM3 
890 CC( I ) = CC ( I ) + (AA( I ,K)*CC(K) ) 

CCl I )=C( I ) — CC ( I ) 
860 CCl I )=CC t I ) /AA( 1,1) 

A(Nt J)-CC(N) 

LIM4=N-1 

CO 910 I=i,LIM4 

I I=N-I 

LIM5=II+1 

CC 920 K=LIM5,N 
920 A{II,J) = A{II,JJ + {AAUI,K)»A(K,J) ) 
910 A{II,J)=CCIII)-AJII,J) 

IFIN-J)110,110,120 
120 C(J)=0. 

CU+l) = l. 

CO 130 I=1,N 
130 CC(I)=0. 

GO TC 840 
110 PUNCH 5 



CC 160 1=1,40 
160 NC(I)=I 

CO 170 J=1,N,5 

PUNCH 6,NCU>iNGU+l)»NCU + 2) ,NCU + 3),N0(J + 4) 
CO 17C 1=1, N 

PUNCH 7,NG( I ), At I , J ) , A( I , J+l) , M I , J + 2) , A( I , J + 3) , A ! 1 , J + 4 ) 
17C CONTINUE 
210 IF( INXT)296,296,99 
296 PAUSE 

CUM=EXITF( 1. ) 

900 PRINT 3 
PRINT 1C 
PUNCH 1C 

IF(LSLN)210, 210,901 

901 CO 9C2 1=1, N, 6 

REAC 2,Ct I ),Ct I+1),C( 1+2) ,C( I+3),C< 1+4) ,C( 1+5) 
9C2 CONTINUE 
KNO=KNO-1 

IF(KNC)210,210,901 

1 FCRfATl.15,15,15,15,15,15) 

2 FORKAT16E10.4) 

3 FORMAT (55H 

4 FCRCATI//43HVALUE OF DETERMINANT OF COEFFICIENT MATR I X = E 1 1. 5 ) 

5 FORK AT (//31H INVERSE OF COEFFICIENT MATRIX — ) 
"6 FCRPAT (//9XI3, 1CXI3, 10X13, ICXI3, 10X13) 

7 FORI" AT{ I3,2XEll.5,2XEll.5,2XE11.5,2XEll.5,2XE11.5> 

8 FCRCAT I //41H SOLUTION VECTOR ACTUAL C-CALC. C) 

9 FCRf AT(2HX( I3.4HJ Ell. 5, 7 XE 11. 5) 

10 FCRCATI/24HINPUT MATRIX IS SINGULAR) 

11 FORMAT ( 17H*CRCUT REDUCTION*//) 
END 
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•CROUT REDUCTION* 



SAMPLE PROBLEM (FOURTH ORDER MATRIX) 



VALUE OF DETERMINANT OF COEFFICIENT MATR I X= .54C00E+02 



SOLUTION VECTOR ACTUAL C-CALC. 

XI 1) .10000E+01 .0C0C0E-99 

X( 2) -.lOOOOe+Ol .00000E-99 

X( 3) .20000E+01 .0CCC0E-99 

X( 4) -.2C000E+01 .O0O00E-99 



INVERSE OF COEFFICIENT MATRIX — 



1. 2 3 4 5 

1 -.33333E-00 .55555E-00 -.14314E-0C .37037E-01 .00000E-99 

2 .83333E-00 -.55555E-00 .14814E-00 -.37037E-01 .OOCOOE-99 

3 .66666E-00 -.44444E-00 -.14314E-00 .37037E-01 .00000E-99 
A .16666E-00 -.11111E-00 -.37037E-01 .25925E-00 .000OOE-99 



